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Abstract 

Properties of the pion are studied at finite temperature with the help of 
the PCAC and the QCD sum rule. The pionic mode is treated consistently 
with the thermal pions which consist of the ground state at finite tempera- 
ture. The temperature dependence of the pion decay constant is estimated 
for niq = 0. No consistent solution is found above the critical temperature 
Tc = \/2fJ^^. The QCD sum rule shows that the finite quark mass makes 
the critical temperature move up to a temperature where the dilute pion gas 
approximation is no more valid. It also suggests that the pion decay constant 
and the continuum threshold decrease by ~ 15% and ~ 10%, respectively, 
and that the pion mass slightly (~ 3%) increases at T = 160MeV, where the 
pion gas approximation seems to break down. 
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1 Introduction 



Behaviors of the hadron properties at finite temperature have been studied exten- 
sively in the context of the phase transition of QCD. In particular, the temperature 
dependencies of hadron masses, coupling constants and so on have been calculated 
in various approaches One of them is the QCD sum rule generalized for finite 
temperature by the authors of ref.|0. 

The QCD sum rule is based on the operator product expansion (OPE) of a 
current correlator and the hadron duality [Q. At first |]^, the Matsubara Green 
function was applied for the massless quark propagator in the OPE of current 
correlators at finite temperature. It was pointed out, however, that the use of 
the finite temperature Green function is not fully consistent for the calculation at 
temperature below the critical temperature, Tc P, |3, §. In the hadronic phase quarks 
have discrete spectrum due to confinement, so that the use of the finite temperature 
Green function requires the full range of their interactions. 

On the other hand, OPE is supposed to achieve a factorization of the soft scale 
and the hard scale. The soft scale dynamics is contained in the matrix elements of 
local operators, while the hard scale dynamics is taken into account in the Wilson 
coefficients. Therefore it is natural to include the temperature dependence only 
in the matrix elements of operators as far as T ~ ^qcd Qi where Q is the 
Euclidean momentum of currents. We conclude that construction of a finite tem- 
perature ground state is essential for the QCD sum rules at finite, but not high, 
temperature. 

In this paper we investigate the pion at finite temperature, since the pion is the 
dominant particle consisting of the finite temperature ground state. At temperature, 
T < Tc, the particles constructing the heat bath (finite temperature ground state) 
are hadrons, especially the pion, since the heavier states are suppressed by the 
Boltzman factor e"™/^ at low temperature. 

We employ the thermo field dynamics of pions in describing the ground state |^ 
and show that the PCAC is valid even at finite temperature. The finite temperature 
ground state |T) gives clear understanding of the temperature dependence of the 
state and naturally avoids a problem of the analyticity at finite temperature 
which was considered to unjustify the QCD sum rule approach. 

The outline of this paper is as follows. In Sect. 2, we define notations and con- 
struct the QCD ground state at finite temperature as a pion gas in the frame work 
of thermo field dynamics. In Sect. 3, we consider the PCAC in the ground state at 
finite temperature constructed in Sect. 2. In Sect. 4, the temperature dependence of 
the pion decay constant is investigated in the chiral limit. There we use the di- 
lute pion gas approximation to estimate the finite temperature matrix elements. In 
Sect. 5, the effects of the finite quark mass to the critical temperature in Sect. 4 and 
other physical parameters of pion are studied by the QCD sum rule technique. In 
Sect. 6, the conclusion and discussion are given. 
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2 QCD ground state at finite temperature 



We define the creation and annihilation operators of the pion, a°^(p) and a"(p), 
which create or annihilate a pion with isospin a and the three dimensional momen- 
tum p. These operators satisfy the following relations: 

= 6\q-p)6^' (1) 

a{p)\n) = y/n\n — 1) 

a^{p)\n-l) = v^|n) (2) 

where \n) is a state that contains n pions with momentum p. The free pion field is 
defined by 

[a"(fc)e-^^" + at"(fc)e^^"] (3) 

■ /(27r)32A;o 



where k = {ko, k) and ko = y k"^ + m^. 

We express the ground state(vacuum) at finite temperature as |T), which is 
regarded as a gas of free pions. Then the free pion annihilation operator cannot 
annihilate |T). Here we follow the thermo field dynamics formulated by Umezawa 
et.al. in ref.|p. They introduce the Bogoliubov transformation and define a new 
operator d;"l"(fc) by 

a{k) = cosh6' a{k) + smh6 a\k) (4) 

Eq.(^) is canonical transformation. The mixing angle 6 = d{k,T) depends on the 
momentum k and the temperature T. 

The tilde operator annihilates a hole, or creates a particle in a state occupied 
in the finite temperature vacuum. In eq.(§) the annihilation operator in the true 
vacuum is divided into an annihilation operator of a particle and a creation operator 
of a hole in the finite temperature vacuum. 

Now we define the finite temperature vacuum |T) such that 

a{k)\T) = 0, &{k)\T) = (5) 

for all k. 

We determine 6 by calculating the ground state average of the particle number 
density which must coincide with the Bose-Einstein distribution, i.e. 

nB{k,T) = {T\a\k)a{k)\T)/6\0) = (sinh^)^ (6) 

with nsiuj) = l/{e^^ — 1) and (3 = 1/T. Then we get 



cosh6' = (1 sinh6' 



n 



1/2 



As we have argued in Introduction, it is regarded that the quark degree of free- 
dom does not explicitly depend on temperature in the confined hadron phase. This 
assumption means that the temperature effect is only in the non-perturbative contri- 
butions of QCD such as the quark condensate. This is of course valid only when the 
temperature is below the deconfinement phase transition. We assume furthermore 
that |T) consists only of pions, since heavier degrees of freedom are suppressed by 
e~™^ at low temperature. 
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3 PCAC at finite temperature 

In order to study properties of the pion at finite temperature, we consider the 
following currents, 

Jt"" = qi.l^t'^q (7) 

Jt = q^l^t'l (8) 



where = ^ is the isospin operator. In the chiral limit the axial-vector current is 



conserved, d^J^"" = Iniqqi'y^t'^q — > 0. In following we omit the isospin indices a for 



simplicity. 

It is important to note that the time-ordered correlator is not analytic in the 
plane at finite temperature |]10[. Instead, the pionic mode can be studied in the 
retarded correlation function, 

n(g) = Jd'xe''^^td,{T\R[J^{x),MO)]\T) 
= q,Jd'xe'^^T\R[J^{x),MO)]\T) 

= qJd^xe^'^^Y.dii) 

'e-'^P--P^>{T\J^\nm){nm\J^\T) - e'^^-'^^'^'' {T\J^\nm) {nm\J^\T) 

(9) 

Here we insert the complete set 1 = X^nm \nfh){nfh\ as the intermediate state. \nm) 
denotes a n-particle m-hole state and is the 4-momentum of the n-particle state, 
while — is the 4-momentum of the m-hole state. 

The PCAC equation, ^^'J^ = IrUgJ^, and 9'^(T| J;f |nm) = -i{pn-pn,Y{T\J^\nTh) 
leads to 

{T\jK,\n'm){nm\J^\T) = — (T| J^lrim) (nm| J5IT) 

When we take the rest-frame, i.e. q = and go = the above correlation function 
is given by 

X {2nf 5(a - (E„ - Em))6^{p^ - pj (T| 4\nm) {nm\J,\T) 

n,m 

= i / dap{a)- — —— (10) 

J [00 + ley — (j^ 

On the other hand, n(co') can be evaluated diredtly, giving 

n(o;) = j d'xe'-Hd,{T\R[J^{x),UmT) 

= J d'xe'^H6{t){T\[4{x),M0)]\T) 

= (T|2gg(0)|T) (11) 
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Thus we obtain 

2uj 



dap{a)-—-- = 2{T\qq\T) (12) 

[uj + ley — cr^ 

Eg ■ (IT2D is valid for any u, i.e. even at u; ^ and the RHS does not depend 
on uj. Then LHS also should not depend on uj. Therefore, we conclude that if 
(T|gg(0)|T) 7^ 0, this correlator is saturated by a massless pole, p{a) = AS^a"^). 
Other states do not contribute to this correlator at all in the chiral limit. 

The spectral function p{a) is defined in eq.(p!0|), 

pia) = i2nfJ25i^ - (En - E^))5'ip^- pJ{T\J,^\nm){nm\J,\T) (13) 

m,n 

The finite temperature vacuum satisfies {T\J\T) = for J = J5 or Jq because of 
parity conservation ^. Then the lowest energy contribution to the intermediate states 
comes from \nm) = |10) and |01). Here we remember that the tilde state is a hole 
and has the opposite quantum numbers of a corresponding particle state. Therefore 
the condition p(cr) = AS^a"^) indicates that the two states |10) and |01) represent 
a massless particle and a massless hole, respectively, and are the only states that 
contribute in the chiral limit. These are the pion and its hole at finite temperature. 

One might wonder why there are two Goldstone modes at finite temperature. 
This is caused by the existence of the tilde current which also satisfies the current 
conservation law. 

At zero temperature a massless pole in the correlator requires existence of an 
asymptotic massless field, J^} -^t-^oo —V^f-Kd^^. Although the finite temperature 
vacuum is not Lorentz invariant , we expect that the same relation for is valid 
in the rest frame at finite temperature except that the decay constant may depend 
on temperature, i.e. ^t-*oo —V^fn^o^- Now we find 

Uico) = 2{T\qq\T) 

= 22(27r)3 ((T| Jo^|10)(10| JslT) + (T| Jo^|Oi)(Oi| J5IT)) 

= 2t{27rf[{T\ - v^/J^o'^|10)(10 r " "^1^) + 



2mq 



^ 2(/J)^(m r^ 
2mq 

2mg 

On the other hand, the retarded correlator of the free pion at finite temperature 
is given by 

{T\R[doHx),<l>m \T) 



{{ub + 1) - Ub} 



^ The finite temperature ground state is connected to the true vacuum by |T) — G ^(6')|0) with 
G{e) = exp [/ <Pke[a{k)h{-k) - at(fc)at(_fe)]] . 
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-(-«)^(xo)(K + l)-nB) 




^—ik()t^—ikx 



ko + ie + Jk +ml 



\T) 



(27r)42 (A;o + le) 



(15) 



Clearly the final form does not depend on the temperature and has the same expres- 
sion as that at zero temperature. It should be noted that the sum of the strengths 
of a particle and a hole state is normalized to the zero temperature pion propagator. 
The difference between eqs. (0) and ( plSf ) comes from the temperature dependence 
of U in the PCAC relation, Jq = -V2/j9o$. 

Finally we obtain the finite temperature version of the Gell-Mann-Oaks-Renner 
relation. 



2mq 



iT\qq\T) (16) 



This relation is same as that at zero temperature except that LHS means the sum 
of a hole and a particle contributions. 



4 Pion decay constant, /J, at finite temperature. 

We estimate the quark condensate at finite temperature in the dilute pion gas 
approximation. As far as T is not high, the expectation value of an operator, O, at 
finite temperature, i.e. (T|0|T), is given ^ by 

(r|0|r) ~ (0|O|0) + ^ / d'p{T:\p)\0\T:'^{p))nB{e) (17) 

a=l 

where e = ^Jp^ + m^, and a denotes the isospin index. To estimate (7r|0|7r) we use 
the soft-pion theorem. Since pions which consist of the pion gas are also affected by 
the heat bath, the second term in eq (|17D, i.e., the matrix element, (vr'^(p)|gg|7r"(p)), 
must also be temperature dependent. As we have seen in eq.(^) the free one-pion 
state Ivr) = a^\Q) does not depend on the temperature. We conjecture that this effect 
can be taken into account by changing the reduction formula from a) d^J^ 

to a"!" '^/ifrd^J^- Then we get 



2^2 
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This is easily solved self-consistently under the assumption m^/rn^ ° 
we find the temperature dependence of the pion decay constant, 



T\2 



cT=0 



2^2 



1 and 



(19) 



The T dependence given by ( pi?D is illustrated in Fig.|I|. It is clear that there exists 
no solution of eq.(|18|) above T > Tc = a/2/J=°. /J decreases towards /J'= = ^^^^ 

when the temperature increases from to y2/J=°. This singular behavior comes 
from the consistency condition that the thermal pions are modified themselves in 
the finite temperature ground state. At T = Tc, the pion gas picture seems to break 
down and there is no more consistent solution above Tr. 



5 QCD sum rules for the pion at finite tempera- 
ture 

So far, we have studied the properties of the pion in the chiral limit, i.e., at 
nig = = 0. We study effects of the finite quark mass to the critical temperature 
and the other physical parameters of the pion by the QCD sum rule technique. 
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5.1 spectral function 

We reconsider the spectral function of the retarded pion correlation function 
defined by 

n^(g) = t J d'xe'^^T\R[J^{x), MO)]\T) (20) 
Similarly to eq.(|^), we find at the rest frame 

n,m 



X 



(21) 

Using a relation of the thermo field dynamics 0, 

{T\J^\nm){nm\J^\T) = e^^^"'^"'\T\J^\mn){mh\J^\T) (22) 
we obtain the following expression, 

= n^=o(c^) 

= ij^^dtj ciae'("-")*(27r)3(l -e-^") 

Y: 5{a - (E„ - Er.))5\p^ - pJ{T\4\nm){nm\MT) 



X 



dapAa) ^- (23) 

uj — a + ie 

where the spectral function is defined by 

pi(a) = (2vr)3(l - e-^"^) 

xY,5{a- (K - E^))5\p:^ - pJ{T\Ji\nm){nm\MT) (24) 

From eq.(P^ and the relation 

{T\J^\nm){nm\jQ\T) = — (T| J(f |?2m) (nm| J5IT) 
we obtain the following symmetry relation for the spectral function pi(cr), 

Pi((T) = pi(-a) (25) 
Then we can rewrite the correlation function. 



LuHiuj'^) = — da2ujpi(a 
Jo 



' * ' ' (cu + ie)^ — 0"^ 



a'^ — io'^ — ie 



^r^^^Ei^ (26) 



Here the last equality is given for the delta function part in the spectral function. 
Eg. ([26|) corresponds to eg. (p^) with pi(cr) = p(cr) + p(— a). From eg.(0), we ex- 
pect that the strength of the pion pole in this spectral function, pi(cr), is given by 

/ rT\1 { T\'2 

2m • Then we may assume the spectral function, pi((T^), in the following form 
for a finite guark mass. 



2m, 



2pi(a^) = ^ 6(a^ - (miy) + Perturbation ■ Oia' - Sq) (27) 



The Perturbation term is estimated by the bare-loop calculation. This form is same 
as that at the zero temperature except that the pion pole consists of a particle state 
and a hole state. 

From our derivation of the spectral function, it is clear that the correlation 
function, n(u;^), has an analytic property in the uJ^ plane. Therefore we use the 
QCD sum rule technigues for this correlation function in the usual manner. 

Here we write down the imaginary part of correlation function, H(u;^). 

-lmH(s) = ^■^"^ ^ 5{s - (mlY) + Perturbation ■ 9{s - sq) (28) 

5.2 The construction of the sum rule. 

Next we construct the theoretical side of the QCD sum rule. The scalar part of 
this correlator can be easily calculated using the Ward identity and the pseudo-scalar 



correlator [|TT 

—OJ 



(2m,){ ^ ^(17 _ ^^^(_^2/^2))| ^ 



4:7ras 
2 _ 



(29) 



At finite temperature the OPE side contains operators carrying spin We 
do not consider operators which have Lorentz indices except for the operators of 
dimension 3, since those operators are proportional to higher powers of temperature 
and also are multiplied by the small guark mass. 

The only possible dimension 3 non-scalar operator is S.T.{d'~^^'^^d + u^^'^^u), 
where S. T. means the symmetric and traceless tensor. This happens to be zero in 
the isospin SU{2) symmetry limit, and therefore no spin operator contributes in this 
case. 

Same as the phenomenological side, we only consider the same theoretical side 
with the zero temperature sum rule except for the temperature dependent conden- 
sate. 
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As is mentioned previously, as far as T is not so high and the thermal pion gas 
is dilute, (T|0(/i^)|T) is approximated as 



{T\0{^^')\T) ^ (0|O(/i2)|0) + E / dMr^^{p)\0{^^)\1x'^{p))ns{e/n 



a=l 



We estimate the temperature dependence of the condensate in the soft pion 
limit [^]. We write down only the results. 

(gg)T = (gg)o(i-^^i?i(^)) 

{G-G)t = {G-G)o 
1 28 ST"^ TTiF 

Here the function Bi{z) is defined by 

6 roc 



Bi{z) = - r dy^y^-z^{ey -r 

TT J z 



-1 



Equating the theoretical side and the phenomenological side of the correlation 
function with the help of the dispersion relation, and making the Borel transforma- 
tion, we obtain the sum rule for the pion at finite temperature: 



2(/J)2(m^)2exp-(-^)'/^^' ^ . 3 ,^ ,17 



2™,|^(l + (^ + 2o^,)^)(l-exp(-4/M=)) 



2 

■]t^(^^)o(l-^5i(^))- (30) 



5.3 The estimation of temperature dependence. 

In its early stages, the QCD sum rule for the pion was known to have a difficulty 
mainly because the coupling of the pion pole is not dominant. For the pseudo-scalar 
correlator, the excited pion states and the instanton effects become important, while 
for the axial-vector correlator, the axial-vector meson Oi will disturb the sum rule. 

These difficulties, however, can be circumvented by adapting the pseudoscalor- 
axialvector off-diagonal correlator since it is almost saturated by the pion. We, 
however, should not take a logarithmic derivative of the lowest moment sum rule to 
eliminate the coupling strength since the logarithmic derivative strongly enhances 
the contribution from the higher resonances. Therefore we evaluate the values of 
physical parameters by fitting the LHS and RHS of the sum rule. 

The fitting analysis of the sum rule is performed in an interval of M, VL{Mmin < 
M < Mmax), the lower and upper limits of which are determined in the following 
argument. 
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According to the Ward identity, the deviation from the chiral hmit is given by the 
pseudo-scalar correlator. The convergence of the OPE of the off-diagonal correlator 
depends on one of the pseudo-scalar correlator. Then we set the lower limit to ensure 
that the power correction of the pseudo-scalar correlator is smaller than 5% of the 
perturbative contribution. This condition is also expected to suppress the instanton 
contribution in the pseudo-scalar correlator. The upper limit is set to suppress the 
continuum contribution to less than 20% of the total in the pseudo-scalar correlator. 

Because Q generally depends on temperature, we determine Q at each tempera- 
ture. At zero temperature the above conditions give at Sq ~ 2GeV^, 

n{0.99GeV <M < l.lSGeV) 

We determine and A(T) = 2{f^f{mlf by matching the RHS, R{sl,T,M), 
and the LHS, L{m^, A(T), M) of the sum rule in the interval Q. Let us rewrite the 
sum rule in the form 

A(T) = F(m^, 4, T, M) = (m, + m,)/2(4, T, M)M^ exp((m^) VM^) 

The coupling parameter A(T) should not depend on M. Therefore we choose 
and Sq so that F is least dependent on M in ^2 and A(T) is determined as the mean 
value F{rn^, s^, T) of the function Flrn^, s^, T, M) in Q. Then we define 

6 = Max[\F{ml, M, T) - F(m^, T)\/F{m^, T)] 

which represents the variation of F{m^, s^, T, M) in Vt. At a given temperature we 
search the best fit parameters which minimize 5. 

Because the condensate depends on the and /J at finite temperature, we 
have to determine and /J self-consistently, i.e., until the difference between the 
input m^, /J and the output m^, /J becomes less than 1%. 

The results are shown in Figures 2-5. We find the similar behavior of /J as 
that in the chiral limit. The consistent solution of /J disappears above the critical 
temperature ~ ISOMeV. At Te, = UQMeV, sq = 1.93 and = 70.9MeV. 
The critical value of the pion decay constant is almost the same as that in the chiral 
limit. It is, however, generally suggested that the pions can no more dominate 



the finite temperature ground state at T > IGOMeV [0. Therefore we trust our 
results only at T < IQOMeV. In this region the pion mass is almost stable, the decay 
constant decreases by~ 15% and the continuum threshold decreases by ~ 10%. 

In this analysis the quark mass and quark condensate are determined by the 
pion sum rules at zero temperature and the 4-quark condensate is evaluated in the 
rho meson sum rule. These values are m„ -|- = 15MeV, {qq)o = — (220MeV")^, 
{qq)l = 7.3 • IQ-^GeV^, and the gluon condensate {^G ■ G)o = {360MeV)\ 

6 Conclusion and Discussion 

The QCD sum rule is applied to the analysis of the pion properties at finite 
temperature. We find that the pion decay constant decreases, while the pion mass 
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Figure 2: The temperature dependence of the pion mass. 
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Figure 3: The temperature dependence of the continuum threshold Sq. 
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Figure 4: The temperature dependence of pion decay constant /J. 
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Figure 5: The ratios m^/m^ °(sohd hne), Sq /sq ° (dashed) and the fj/fj ° (dot- 
ted). 
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increases as the temperature increases. These results quahtatively coincide with 
those in other approaches, such as the NJL model ||13|. The increase of the pion 
mass is very small, i.e., ~ 3% at T = IGOMeV, which is again consistent with the 



result in the chiral perturbation theory [14]. 



It is found that the consistency of the pion properties with those of the thermal 
pions suggests a critical behavior of the pion decay constant. In the chiral limit the 
critical temperature, is y2/J=° ~ 133MeV for /J=° = 9AMeV. Above T^, there 
exists no consistent solution of /J. The higher power corrections of T for {qq)T can 



be expressed in the chiral limit as follows [12 



mr = moil - ^ - ^ - ^^^4 + ^^'^ 

One sees that they tend to accelerate the breakdown. Furthermore the contribution 
of the higher resonances is expected to suppress the quark condensate. Therefore 
these corrections to our analysis are expected to enhance the critical behavior of the 
pion decay constant. 

Our results of the "phase transition" seems to be the first order phase transi- 
tion. The QCD phase transition with two massless quarks is expected to be the 
second order [0. We conjecture that the discrepancy comes from the breaking of 
hadronic description of the finite temperature vacuum. At the temperature near 
Tc(= V2/J=°), the quark de grees of freedom are no more frozen, and the self-energy 
of quarks will break the PCAC relation. Such an effect might make the "phase 
transition" smooth. 

The QCD sum rule with a finite quark mass indicates the same critical behavior 
at T = ISOMeV. The difference of the critical temperature from the chiral limit can 
be understood from the consideration that the small difference of the quark mass 
leads to a large (compared with the quark mass) scale difference in the hadronic 
world, such as the pion mass. This, however, may not justify the hadronic descrip- 
tion of the finite temperature vacuum at T = ISOMeV, since the higher resonance 
contribution is expected to dominate the number density at T > IGOMeV. 

It is worthwhile to mention about the relation between Sq and {qq)T- The pion 
mass is approximately calculated by taking a logarithmic derivative of eq.(^). If 
we omit a few negligible terms and take M ~ IGeV then we get 



= 



-2{qq)T 

When the pion mass increases at finite temperature this equation gives the following 
constraint: 
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T 



-2{qq)Timi=y < 2m,— [1 - (1 + si )e-^o 

If we regard the continuum threshold sq as an order parameter for the decon- 
finement phase transition [|l^, then the above constraint indicates that the critical 



temperature T^, defined by Sg" = 0, cannot be lower than the chiral transition 
temperature Tch with (gg)^^^^^ = 0. If the pion mass does not depend on the 
temperature, these two temperatures are equal in our estimation. 
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